Unconventional Superconductivity from Local Spin Fluctuations in the Kondo Lattice 
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The explanation of heavy-fermion superconductivity is a long-standing challenge to theory. It is commonly 
thought to be connected to non-local fluctuations of either spin or charge degrees of freedom and therefore of 
unconventional type. Here we present results for the Kondo-lattice model, a paradigmatic model to describe 
heavy-fermion compounds, obtained from dynamical mean-field theory which captures local correlation effects 
only. Unexpectedly, we find robust s-wave superconductivity in the heavy-fermion state. We argue that this 
novel type of pairing is tightly connected to the formation of heavy quasiparticle bands and the presence of 
strong local spin fluctuations. 

PACS numbers: 71.27.+a,74.20.-z,74.20.Mn 



Heavy-fermion (HF) 4f and 5f intermetallic compounds 
constitute a paradigm for strong electronic correlations. Their 
low-temperature behavior is affected by f-shell local moments 
subject to antiferromagnetic (AF) exchange coupling to the 
conduction electrons, resulting in Fermi-liquid (FL) phases 
with strongly renormalized Landau parameters, most notably 
huge effective masses [1-4]. Moreover, HF materials often 
display symmetry-breaking phases which occur either within 
the heavy FL [1] or compete with it [5]. While magnetic order 
in systems containing unscreened moments appears natural, 
HF superconductivity is conceptually non-trivial and indeed 
came as an unexpected discovery more than three decades ago 
[6]. By now, a wide variety of f-electron superconductors are 
known [7, 8], many of them confirmed to be unconventional 
[1,2, 9]. 

Superconducting (SC) transitions in HF compounds are of- 
ten assumed to be driven by non-local fluctuations of the f- 
shell spin degrees of freedom. This idea finds support in the 
close connection between HF superconductivity and magnetic 
quantum phase transitions (QPTs) where such fluctuations are 
strong [3, 5, 10, 11]. Alternatively, pairing mediated by fluc- 
tuations in the charge channel (i.e., valence fluctuations) has 
also been discussed [12]. Given that the basic theoretical mod- 
els for HF materials, the periodic Anderson model (PAM) and 
Kondo-lattice model (KLM), constitute complicated interact- 
ing many-body problems which cannot be solved exactly, the- 
oretical descriptions of HF superconductivity often employ ei- 
ther simple static mean-field theories or effective models of 
fermions coupled to spin or charge fluctuations. 

A rather successful approach to study the microscopic 
properties of correlated-electron lattice models beyond static 
mean-field or effective descriptions is the dynamical mean- 
field theory (DMFT) together with its cluster extensions 
[13, 14]. The lattice problem is mapped onto a self-consistent 
quantum impurity model at the expense of losing informa- 
tion on non-local correlation effects beyond the spatial size of 
the impurity cluster. Therefore, it is commonly assumed that 
a proper description of HF superconductivity within DMFT- 



based approaches requires either large enough clusters or the 
inclusion of a bath in the two-particle channel (which explic- 
itly models a bosonic "glue" for superconductivity). In partic- 
ular, within the conventional (single-site) DMFT only s-wave 
superconductivity occurs [15] and a relation to magnetic fluc- 
tuations appears highly unlikely. 

In this Letter we report on the unexpected observation of a 
stable SC solution to the DMFT equations for the KLM with- 
out any external glue. Although the pairing symmetry is s- 
wave, the SC state is highly unconventional: Pairing is driven 
by local spin fluctuations; it comes with a strong frequency 
dependence of the gap function and requires the formation of 
HF bands as a prerequisite. 

We note that a hint of the possible occurrence of local su- 
perconductivity was found in an earlier DMFT study to the 
PAM [16] which, however, did not analyze the SC phase but 
only normal-state instabilities. Static mean-field descriptions 
of the KLM or the PAM can also yield solutions with local 
pairing [17-19], but it is difficult to assess their validity, as 
fluctuations beyond mean field may destroy pairing. 

Model. Within the KLM it is assumed that the f-states are 
sufficiently localized and can be described by a (pseudo-)spin 
degree of freedom which couples to the conduction (c) elec- 
trons via an exchange interaction. We use the simplest version 
of the model, i.e., a nearest-neighbor tight-binding conduction 
band with spin degeneracy only and a S = 1/2 spin located 
at each lattice site. The Hamiltonian reads 

H = "* E £ l £ ja + ^ES r &T Q/J C i/} . (1) 
(i,j),<7 i-afj 

Here cfj denote annihilation (creation) operators of conduc- 
tion electrons with spin a at site i, and (., .) denotes nearest- 
neighbors. Si is the operator for the localized spin, t the vec- 
tor of Pauli matrices, and the interaction between the conduc- 
tion electrons and the localized spin is modeled as an isotropic 
exchange coupling with J > 0. For simplicity, we will con- 
sider nearest-neighbor hopping t on the infinite-dimensional 
Bethe lattice, leading to a semi-circular density of states with 
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bandwidth W. We have checked that using different lattice 
types (for example hypercubic or square lattice tight-binding) 
does not change the results qualitatively. 

Methods. Within standard DMFT, the conduction- 
electron self-energy is approximated as local in space, 
E(k, ui) — > E(u;). Then, the KLM maps onto an effective 
single-impurity Kondo model (SIKM) [4] which needs to be 
solved within a self-consistency loop. The effective SIKM 
couples to the dynamical mean field, i.e., a bath which is char- 
acterized by the Green's function Q^ 1 (ui) = (uj) + E(w), 
where G\ oc denotes the local effective impurity Green's func- 
tion, and E the self-energy. There are several computational 
tools to treat the SIKM, for instance the continuous-time 
Monte Carlo [20, 21] or Wilson's numerical renormalization 
group (NRG) [22]. In this work we utilize the NRG, its 
advantages being that one can access arbitrarily small energy 
scales, calculate spectra directly on the real-frequency axis, 
and work at both T = and T > 0. We work with the dis- 
cretization parameter A = 2.0, keep N st = 1000, . . . , 2000 
states, and perform z-averaging with N z = 2 [22]. 

To allow for solutions with SC order, we generalize the 
DMFT equations and the impurity solver to a Nambu formula- 
tion with 2x2 matrix propagators [22, 23]. Together with the 
use of DMFT, this constrains our calculations to spin-singlet 
even-frequency s-wave superconductivity. We recall that the 
DMFT treatment of superconductivity is non-perturbative and 
thus goes beyond the standard Eliashberg theory [24]: it does 
not rely on any assumption about a separation of energy scales 
for the fermions and the bosonic glue responsible for the for- 
mation of superconductivity. At present, we restrict the calcu- 
lations to SC order only, suppressing possible magnetic order. 
Our result below show that strong pairing occurs in a regime 
without magnetic order, justifying this neglect. 

First results applying the DMFT+NRG method to super- 
conductivity in the Kondo lattice were published in Ref. [25], 
where an additional negative-?/ Hubbard interaction term for 
the c electrons was used to induce pairing which then com- 
peted with the Kondo effect. In the present work, we consider 
instead the KLM without explicit attractive interaction in the 
Hamiltonian; a SC solution had not been expected in this case. 

Results: Superconductivity at T — 0. Our numerical solu- 
tion of the DMFT equations yielded, for a range of model pa- 
rameters, stable SC solutions whose properties we discuss in 
the following. The conduction-band density of states (DOS) 
of both the paramagnetic normal (N) and SC solutions of the 
KLM for fixed conduction band filling n = 0.9 and different 
J are shown in Fig. 1, left panels. The only feature of the N 
DOS is a hybridization pseudo-gap above the Fermi energy, 
signaling the formation of heavy quasiparticles, which can be 
rationalized within the picture of hybridized c and / bands. 
The N solution is unstable against superconductivity, where 
the DOS exhibits two additional features: (i) A true gap A sc 
with well-developed van-Hove singularities is present around 
the Fermi energy; as a function of J, it first increases up to 
J/W — 0.5, and then slowly decreases, (ii) In addition to 
the SC coherence peaks, there are side resonances at positions 
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FIG. 1: (color online) Left panel: N DOS (dashed/red lines) and SC 
DOS (full lines) for n = 0.9 and various J. Right panel: Real part 
of the corresponding gap functions. 



which roughly scale with J. These structures are sharp for 
smaller J, but become increasingly washed out for larger J. 
These features are likely related to local spin fluctuations sta- 
bilizing the pairing, as discussed in the supplemental material. 

As the appearance of a gap alone is not sufficient to iden- 
tify the solution as a SC, one needs to look at the anomalous 
parts of the Nambu Green's function respectively the anoma- 
lous part of the self-energy. From it, a SC gap function can be 
defined, as in standard Eliashberg analysis, via 
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where E Q (w) denote the components of the electronic self- 
energy expanded into Pauli matrices, E = E q Tq, {a — 
0, 1, 2, 3), in Nambu space. The resulting real parts ReA(cj) 
are shown in the right panels of Fig. 1. As expected for even- 
frequency pairing, ReA(w) is symmetric. It shows a strong 
frequency dependence, with sharp features shifting to larger 
energies and broadening with increasing J. These structures 
are linked to the side resonances in the DOS: the zeroes of 
ReA coincide with the local minimum of the resonances in 
the DOS. The uj — limit, ReA(0), provides an estimate of 
the gap seen in the DOS; it exhibits the same non-monotonic 
behavior with J as noted above for the gap in the DOS. 

To characterize the evolution of the superconductivity 
across the phase diagram, we plot in Fig. 2 the anomalous 
expectation value $ = (cj+Cji ), as a function of ,/ and n. Su- 
perconductivity is found to be stable over wide regions of the 
phase diagram for J/W > 0.1. For a fixed J/W = 0.2, a 
finite $ is found between 0.45 < n < 1. For larger J/W the 
SC region extends to even lower fillings. A maximum of $ 
appears around J/W = 0.3 and n — 0.9; the side-resonances 
are also the most pronounced for these parameters. 

In Fig. 3 we display the evolution of the order parameter 
$ along two cuts along the phase diagram indicated by the 
white dashed lines in Fig. 2: for fixed n = 0.9 (Fig. 3a) and 
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FIG. 2: (color online) Anomalous expectation value $ as a func- 
tion of J/W between quarter and half filling. The white dashed 
lines indicate the cuts along a fixed J resp. n shown in Fig. 3 (a), (b). 



for fixed J = O^iy (Fig. 3b). The analysis for weak Kondo 
coupling, J/W < 0.1, is difficult as the signatures of SC be- 
come very weak and hard to distinguish from numerical noise. 
Thus we cannot decide whether the SC solution ceases to ex- 
ist for small J, or whether it survives down to J — > with 
an (exponentially) small pairing scale. (The latter would be 
expected in the weak-coupling limit of certain mean-field the- 
ories, see supplemental material [19].) For J > W/2, on the 
other hand, we observe a decay consistent with $(</) oc 1/J. 
We will comment on this behavior further down. 

We note that even for the half-filled case n = 1, where 
the normal-state solution is a Kondo insulator, the SC solu- 
tion exists in a small region 0.05 < J/W < 0.15. The 
corresponding $ is, however, roughly an order of magnitude 
smaller compared to n < 1. 

Normal-state Fermi-liquid scale. In the normal state, the 
KLM realizes a heavy FL at low temperatures for n ^ 1, with 
a FL (coherence) scale To (typically smaller than the single- 
impurity Kondo temperature Tk)- In a local self-energy ap- 
proximation, To can be efficiently extracted from the quasi- 
particle weight 



dReE(w) 



du 
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via T = WZ, where S(o;) = S (w) + S 3 (w). The evolution 
of To with J and n is also depicted in Fig. 3(a) and (b); we re- 
call that for J — > the scale To depends exponentially on both 
J and the bare c DOS, T oc yj J/W exp (-a(n) • W/J) 
with a weakly n -dependent coefficient a{n) [26, 27]. For 
large J > W/2 the dependence of T on J significantly de- 
viates form this Kondo form and rather tends to saturate as 
J — >• oo. Finally, for fixed J and varying n we recover the 
known dependency To oc n ■ e c ' n [27]. These different types 
of behavior for T)(n, J) can be seen from the lines superim- 
posed to the data in Fig. 3. 

Apparently, there does not exist a simple connection be- 
tween To and $. For small J Q(n, J) seems to scale with 
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FIG. 3: (color online) $ (circles) and To (squares) as a function of 
J at fixed n = 0.9 (a), or as a function of n at fixed J/W = 0.2 
(b). The full lines represent approximate dependencies of To on J 
respectively n (see text). 



To. However, as noted before, the results for very small $ be- 
come unreliable for numerical reasons, i.e. one cannot readily 
extract a simple relation between $ and T in this limit. For 
large J at n = 0.9, on the other hand, we do not see a direct 
relation between To and $, but find $ oc 1/ J instead. 

Competition with magnetism. Within the DMFT for the 
KLM, one also finds magnetic phases, namely antiferromag- 
netism (AF) close to half filling and ferromagnetism (FM) at 
small filling [21, 28, 29]. Note, however, that these phases 
have limited extent in both J and n, for example at n = 0.9 
we find AF only for J < J c (n = 0.9) » 0.2W [21] and FM 
only for n < n c (J/W = 0.2) sa 0.65 [29]. Thus, the re- 
gion in parameter space where we find a strong superconduct- 
ing phase in Fig. 2 seems to be complementary to the regions 
with magnetic phases. As the boundaries seem to overlap - 
in particular the SC phase at n — 1 lies inside the AF regime 
- it is surely interesting to study the competition respectively 
interplay of AF, FM and SC in detail. This is work in progress. 

Results for T > 0. Figure 4 displays finite-temperature 
results for the DOS: With increasing temperature T the gap 
shrinks and the spectral side resonances are depleted. The 
reduced pair correlations are also reflected in a decrease of $ 
(see inset to Fig. 4). Close to T c the gap is progressively filled, 
and both the hybridization gap and the side resonances move 
towards the Fermi level. Finally, in the normal-state solution 
for T > T c only the hybridization pseudo-gap is visible. 

Estimating T c from the numerical data at finite T is hard, as 
close to the critical temperature the signatures of SC become 
very small and also the convergence of the DMFT rather slow. 
From the data in Fig. 4 we extract T c s» 0.0036W / at n = 
0.9 and J = 0.25W, which is well below T . In Tab. I we 
collect the resulting estimates for T c for fixed n — 0.9 and 
several values of J in the region where according to Fig. 3 
we have optimal conditions for pairing. As a rule we observe 
that always T c < Tq, i.e. the HF state seems to be a necessary 
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FIG. 4: (color online) Evolution of the DOS with increasing tem- 
perature for n = 0.9 and J = 0.25W where T /W = 0.0105 at 
T — > 0. Inset: $ as a function of temperature. Pair correlations are 
fully suppressed for T > T c . 
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0.2 


0.25 


0.3 


0.4 


0.5 


T /W 


0.0050 


0.0105 


0.0165 


0.0285 


0.0387 


$(T = 0) 


0.0160 


0.0163 


0.0174 


0.0153 


0.0140 


ReA(0)/W 


0.0138 


0.0165 


0.0180 


0.0193 


0.0195 


T c /W 


0.0027 


0.0036 


0.0054 


0.0058 


0.0054 


T c /T 


0.54 


0.34 


0.32 


0.20 


0.13 


ReA(0)/T c 


5.169 


4.583 


3.321 


3.305 


3.585 



TABLE I: Quantities characterizing the superconducting solution for 
n — 0.9 as a function of J in the region of optimal conditions for 
pairing. 

ingredient for the appearance of the SC phase. For small J this 
obviously leads to a strong suppression of T c . On the other 
hand, the gap ReA(0) appears to be much less sensitive to J 
respectively T in this regime. An interesting characteristic 
quantity is the ratio between the gap and T c . The results are 
given in the last row of Tab. I. Obviously, the ratio exceeds the 
BCS value A/T c w 1.74 by a sizeable factor between 2 and 3. 
It is interesting to note that such values are actually observed 
in HF superconductors [11], although the interpretation there 
is usually given in terms of a weak-coupling theory for a d- 
wave state. 

Pairing mechanism at strong coupling. In the strong- 
coupling limit, J > W, the pairing mechanism and behavior 
of $( J) can be understood perturbatively, Fig. 5. To this end, 
we consider a conduction-band filling of n = N e i /N s < 1 on 
N s sites. For J/W — ► oo there are N e i Kondo singlets and 
(N s — N e i) uncompensated local moments. The c electrons 
are mobile, such that the uncompensated moments can alter- 
natively be interpreted as spinful c holes with density (1 — n) 
and a hard-core repulsion, forming the Fermi liquid with a co- 
herence scale To oc W ( whereas the impurity Kondo scale 
simply diverges oc J). For large but finite J/W additional 



excitations out of this manifold are allowed. The lowest one 
consists of converting a singlet into a triplet and can be cre- 
ated by c (hole) hopping. Interestingly, the excited triplet may 
decay via a different neighboring hole provided that its spin 
is opposite to the first one. Together, this second-order vir- 
tual process leads to correlated hopping, with an energy gain 
oc W 2 / J, which binds two holes into a singlet state, Fig. 5. 
The pairing is local - it occurs on the site of the virtual triplet 
- but the holes share this site only in the virtual triplet state, 
so that the pairing is strongly retarded. 

It is plausible that this pairing mechanism continues to op- 
erate at smaller J. Importantly, the existence of the virtual 
state, whose energy may now be approximated by the Kondo 
binding energy T , requires Kondo screening to be intact - this 
naturally explains the limitation T c < Tq. The picture also 
makes clear that superconductivity is more favorable close 
to half filling: Here, Kondo screening is done by c electrons 
whereas the paired carriers are c holes (for n < 1). In contrast, 
in the opposite (exhaustion) limit of small n, Kondo screen- 
ing becomes strongly non-local, and it is the same c electrons 
doing the screening that need to be paired, thus the pairing is 
weak. 

Summary. We have identified a novel mechanism for su- 
perconductivity in heavy-fermion materials: Local spin fluc- 
tuations due to the Kondo exchange coupling can act as a re- 
tarded paring interaction and drive s-wave superconductivity 
in the heavy FL. A particularly interesting feature is the ap- 
pearance of structures in the tunneling DOS at scales related 
to the spin fluctuation spectrum, i.e. well separated from the 
coherence peaks at the gap edges. Such structures have been 
observed for example in recent STS measurements on iron 
pnictides [30]. It would be interesting, to extend such STS 
experiments to systematically study HF superconductors. 

For model parameters relevant to HF systems, T c can be 
as large as several Kelvins, a typical T c scale for existing HF 
superconductors. However, given that pairing in our theory 
is s-wave, it can be ruled out for those materials where the 
existence of gap nodes have been established experimentally. 
More generally, it is an interesting question to what extent 
this mechanism can cooperate or will actually compete with 
SC driven by e.g. non-local magnetic fluctuations. To address 
this point numerical studies beyond DMFT will be required. 
Work along these lines is in progress. 
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FIG. 5: (color online) Second-order process responsible for pairing 
in the limit of large J/W. Shown are the c electron and local- 
moment (f) configurations on adjacent lattice sites; the ellipses de- 
note a singlet bond of two electrons. Here, two hoppings effectively 
move a down-spin hole by two sites across an up-spin hole, with a 
triplet intermediate state (bold/red). 
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LOCAL PAIRING IN MEAN-FIELD THEORY 

Here we summarize how local pairing emerges in the 
Kondo lattice model within mean-field theory [1, 2]. One 
starts by representing the f-shell local moments in terms of 
auxiliary fermions: 



ia'a-V Jia 



fi( 



(SI) 



where f ia describes a spinful (but charge-neutral) fermion de- 
struction operator at site i. This representation requires the 
constraint Y^, a fi a fic = 1 an d displays a local U(l) gauge 
invariance. 

The standard decoupling of the Kondo interaction in the 
KLM [3] is done via a bosonic field bi in the particle-hole 
channel, resulting in an interaction term of the form 



E 



{bi4 a fia + h.c.) . 



(S2) 



Bose condensation of the bi endows the f particle with a phys- 
ical charge and leads to the standard picture of two hybridized 
electronic bands in the heavy Fermi liquid. Alternatively, the 
interaction can be decoupled in the particle-particle channel, 



n- b = Y,{hc\jt-a + h.c). 



(S3) 



Importantly, condensation of the bi alone yields the same 
metallic HF state, due to the local U(l) gauge invariance. 
Superconductivity is generated by a simultaneous condensa- 
tion of bi and bi which induces anomalous expectation values 
among the physical c electrons, (q-j-c^) ^ 0, and this breaks 
the physical U(l) electromagnetic symmetry [4]. The result- 
ing state displays local s-wave pairing with a full gap. 

This mean-field description of Kondo-lattice superconduc- 
tivity via two decoupling channels of the Kondo interaction 
has, however, obvious deficiencies. In particular, in the limit 
J/W — > oo it yields strong pairing on the scale J; this re- 
sult is incorrect as the Kondo interaction only causes spin sin- 
glet formation on this scale (but no Cooper pairing). In the 
opposite limit of small J, pairing in this mean-field theory 
becomes exponentially weak, with T c scaling as the single- 
impurity Kondo temperature T K . 

As an aside, we note that if one starts from an Anderson 
(instead of Kondo) lattice Hamiltonian, the / particles always 



carry a physical charge, such that only a mean-field decou- 
pling in the particle-particle channel is required (in addition 
to the c-f hybridization) to induce superconductivity. Then, 
the resulting pairing may be understood as local c-f pairing 
[1,2]. 



LOCAL SPIN-FLUCTATION THEORY FOR PAIRING 

Let us try to give an argument why for the KLM in DMFT 
the spin fluctuations can be a driving glue for superconductiv- 
ity. To this end we use the approximate form of the pairing 
equations as given in [5], and employ Hewson's renormalized 
perturbation theory [6], i.e. we assume that we have a local 
Fermi liquid with correspondingly renormalized parameters. 
Note that this also means that the Kondo spin is now a dy- 
namic quantity with a fluctuation spectrum living on the order 
of the Kondo scale. We do not use the usual approximation in 
the Eliashberg equations that all fermionic energies can be put 
to the Fermi energy, because all relevant energy scales are of 
the order of T , including the spin fluctuation scale. 

We are working within DMFT, i.e. all momentum depen- 
dencies can be ignored and the k sums for the Green's func- 
tions be performed using the non-interacting DOS. The equa- 
tion for the gap can then be written as 

oo 

£ (12) (c^0+) = J dJ K{uj,uj') -hmG (12) (Lo' + iQ+) 



where the kernel is given by 

oo 

~9 2 x'm 



tanh — + coth — 



For the integral the principal value has to be taken and x's( w ) 
is the imaginary part of the dynamical susceptibility of the 
Kondo spin. The constant g 2 collects the (unknown) coupling 
matrix element and some constant numerical factors. 

At this point one usually assumes that u, oj' <C [i and puts 
w « oj' « 0, which then leads to 



~2 / 

■^-tanh- X (0) 



As x(0) > 0, this kernel is repulsive and hence will not gen- 
erate s-wave even-frequency superconductivity. As already 



2 



mentioned, we here cannot safely assume u>, u>' <C /i, since 
for our Fermi liquid — > fi ^ Tq. We therefore must try to 
analyze the equation somewhat more carefully. To this end 
we nevertheless will assume that we are deep enough in the 
Fermi liquid phase so that temperature dependencies of the 
Fermi liquid parameters are negligible, and that - at least for a 
determination of T c - we can approximate T,^ 12 ^(uj + i0 + ) « 
—A 9( J/2 — oS) E HL The estimate for the cut-off is taken 
from the evolution of the gap function in Fig. 1. Since the gap 
is frequency independent, we can at least put ui = and - 
ignoring the term with coth [5] - find a kernel 



Putting everything together into the gap equation we arrive 



at 



.7/2 

1 = Np J du>K(0,u>) 
o 



1 



1 



u - V lo + V 



(S4) 



as equation to determine T c , where the integral has to be eval- 
uated as principal value integral. 

The kernel contains x's( u )> which is shown in Fig. SI. 
Although the kernel is negative for small oj and would be con- 



= — !~tanh — xs(-w') 
kg 2 , , , ,s 

= -^"tanh — 

As final input we need the anomalous Green's function, 
which is given by 



£( 12 )( 2 



(z + er-E( 22 )(^))(z- 



S( 11 )(z))-S( 12 )(z) 2 



using the independence of the self-energy on momentum in 
DMFT We replace YS 12 \z) by our approximation and also 
set it to zero in the denominator, as we are interested in finding 
T c . Furthermore, the diagonal self-energies are to be replaced 
by the expressions for the local Fermi liquid, i.e. 
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FIG. SI: Xs(w) for J = Q.2W and n = 0.9. The dashed line marks 
the zero. In the inset the result for the integral on the right hand side 
of Eq. (S4) as a function of V, evaluated for T = 0, is shown. 



V 2 



Z + Jx 



where the effective chemical potential fi and hybridization V 
are both O(Tn). Using this self-energy expression in the nor- 
mal phase leads to the well-known formation of heavy quasi- 
particle bands [6, 7]. In the following we further assume 
fj, = 0. With these simplifications we get 



G (12) 



-A 



1 



1 



z-\e j :\~V 2 /z z + \e j :\-V 2 /z 



We can now perform the k-sum. To this end we assume a 
featureless DOS of the band states with value Np and use the 
wide -band limit to obtain 



G {12) (z) = iirANp 



1 



V 2 /z 



sidered repulsive in a standard calculation, the contribution 
from the anomalous Green's function within the approxima- 
tion of a local Fermi liquid due to the Kondo effect is also 
negative for energies less than some constant times T , hence 
the total contribution from the region u G [0, To] is positive. 
There is of course the pole, which has to be treated in the prin- 
cipal value sense, leading to a sign change for u> > V. Here it 
becomes important that x' s (w) decays rapidly on a scale To, 
too. 

As very simple approximation let us assume that x'si 10 ) ~ 
Xs(0) 0(^0 — M), where cu is some cut-off for the local spin 
fluctuations. According to Fig. SI, a suitable cut-off would be 
again O(Tq). We restrict the discussion to T = 0, because if 
we cannot find a solution there, we will surely not find it for 
any T > 0. Under these conditions the integral on the right 
hand side of (S4) then evaluates to 
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7T<? 2 XS(O) 
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1 



1 



w - V oj + V 



irg 2 Xs(0) 



In 
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As long as ujo < V"2 V, the logarithm on the right hand side is 
negative, i.e. the total expression positive and hence one has a 
non-trivial solution to equation (S4). 

Let us now use the numerical data for x' s (w) and see, if the 
above rough arguments are compatible with them. As param- 
eter we vary the effective hybridization between some small 
value V <C To to some large value V ^> Tq. The result is 
shown in the inset to Fig. SI. In both limits, the value of the 
integral is negative, hence no pairing is supported. However, 
in a reasonably large regime with V = O(To) one indeed ob- 
tains a positive value and hence the tendency to form Cooper 
pairs. 

One thus needs two absolutely necessary ingredients - 
hybridized bands with a small energy scale, and a spin- 
fluctuation spectrum that decays on the same energy scale. 
Note that this result also predicts that one has to expect a 
strongly reduced T c for a phonon-mediated superconducting 
heavy-fermion state, as here the kernel enters with a negative 
sign from the outset (see e.g. [5]), and hence the low-energy 
part of the integral now rather tends to suppress superconduc- 
tivity. 

We note that the above discussion - considering the vari- 
ous uncontrolled approximations - can neither be considered 
a proof of the existence of s-wave superconductivity, nor does 



it provide a reliable formula for estimating T c . For instance, 
we have assumed a frequency-independent gap, obviously in- 
consistent with the full numerical results in Fig. 1 . The discus- 
sion nevertheless gives a qualitative foundation for our results, 
namely that, contrary to common wisdom, spin fluctuations 
can lead to s-wave superconductivity, at least under the very 
special conditions met in heavy-fermion systems. 
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